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Abstract 

A spatially flat Friedmann-Robertson- Walker (FRW) cosmological model with 
generalized Chaplygin gas is studied in the context of scalar-metric formulation of 
cosmology. Schutz's mechanism for the perfect fluid is applied with generalized 
Chaplygin gas and the classical and quantum dynamics for this model is studied. 
It is found that the only surviving matter degree of freedom played the role of 
cosmic time. For the quantum mechanical description it is possible to flnd the wave 
packet which resulted from the linear superposition of the wave functions of the 
Schrodinger-Wheeler-DeWitt(SWD) equation, which is a consequence of the above 
formalism. The wave packets show two distinct dominant peaks and propagate in 
the direction of increasing scale factor. It may happen that our present universe 
originated from one of those peaks. The many-world and ontological interpretation 
of quantum mechanics is applied to investigate about the behaviour of the scale 
factor and the scalar field (considered for this model). In both the cases the scale 
factor avoids singularity and a bouncing non-singular universe is found. 

Keywords: quantum scalar-metric cosmology, Saez-Ballester theory, Chaplygin gas, 
Schutz formalism 

1 Introduction 

In recent years an exotic fluid known as Chaplygin gas has well described the reason behind 
the current observation of cosmic acceleration. It was first proposed in [1]. The generalized 
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Chaplygin gas satisfies an equation of state p — where p is the pressure, A is a positive 
constant and 0<a<l[l, 2]. a = 1 corresponds to standard Chaplygin gas [1]. This 
equation of state describes an universe evolving from a non-relativistic matter dominated 
one to a cosmological constant dominated one [2] . The equation of state for the Chaplygin 
gas actually comes from the string Nambu-Goto action in the light-cone coordinate [3]. 
The idea of this exotic fluid was later used widely [3, 4, 5, 6, 7, 8, 9, 10] mainly to find 
a solution to the coincidence problem in cosmology. The isotropic quantum cosmological 
models with Chaplygin gas have been studied very recently [11, 12]. The generalized 
Chaplygin gas is also studied in a modified gravity approach [13] and the geometrical 
explanation for the generalized Chaplygin gas gives the deceleration parameter and the 
age of the universe which are consistent with observation [14] . 

Here we shall study the quantum dynamics of a spatially flat FRW model with Chap- 
lygin gas and a dimcnsionlcss scalar flcld coupled to the metric. This formulation was 
first classically studied by Saez and Ballester [15] and there a satisfactory answer to the 
missing-matter problem in cosmology is given. Recently the quantum description of this 
Saez-Ballester theory in a generalized form is studied [16, 17, 18]. In our case we take a 
polynomial coupling of the form GA^^™" to the metric. We then apply Schutz's formal- 
ism of a perfect fluid [19, 20], which was extensively used in some recent publications 
[21, 22, 23, 18] to the Chaplygin gas. With the canonical methods [24] we get the su- 
per Hamiltonian. We studied the classical dynamics of the system and found that the 
remaining matter degree of freedom can play the role of time. Canonical quantization 
of the super Hamiltonian gives us the Schr6dinger-Wheeler-DeWitt(SWDW) equation. 
Then we linearly superpose the wave functions of the SWDW equation to construct the 
wave packet. In the next step we recovered the classical solutions using the many- world 
and ontological interpretation of quantum mechanics. 



2 The Quantum Model 

The expression for action in this spatially flat quantum FRW model with generalized 
Chaplygin gas can be written as 

[ d^x^/^ [R - F{(P)g''''<PM + 2 / d'xVhhab i^"" + / d^x^ (1) 
Jm JdM Jm 

where hab is the induced metric over three dimensional spatial hypersurface which is the 
boundary dM of the four dimensional manifold M and X"'' is the extrinsic curvature. Here 
units are so chosen that c = IGttG = ^ is equal to one. Pc is the pressure for generalized 
Chaplygin gas. Here F{(j)) is an arbitrary function of a scalar field which is coupled to the 
metric. The first term(with R only) and the third term were first obtained in [25]. The 
generalized Chaplygin gas satisfies the equation of state 
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Here ^4 is a positive constant and 0<a;<l. a = 1 represents the standard Chaplygin 
gas. In Schutz's formalism [19, 20] the fluid's four velocity can be expressed in terms of 
three potentials e, 9 and S (here we are studying spatially flat FRW model so other 
potentials are absent in this model because of its symmetry), 



(3) 



Here h is the specific enthalpy, S is the specific entropy, e and 6 have no direct physical 
meaning. The four velocity also satisfy the normalization condition 



u'^Ui, — 1. 

The metric for the spatially flat FRW model is 

ds^ = N\t)dt^ - a\t) \dr^ + r\d^^ + sin^ ^dip^] 



(4) 



(5) 



where N{t) is the lapse function and a{t) the scale factor. If we study the thermodynamics 
of a perfect fluid [24] we can see that h and S satisfy these thermodynamic relations 



p = Po(i + n) , h = i + u + p/po , Tds = du + Pd{i/po) 



(6) 



where r and p are temperature and particle number density respectively. We can use 
Schutz's mechanism with these thermodynamic relations to evaluate the third term of the 
action. If we eliminate the surface terms , the first and second term of the action gives 



dt 



So the gravitational Lagrangian density is 

6 1 
= --ah^ + -F{cj>)a^''. 

The conjugate momenta for these a and are then 



(7) 



(8) 



Pa 



and the Hamiltonian for this part is 



nn = N 



"24 a ^ 4F((/))a3' 



(9) 



(10) 



Let us now find the Hamiltonian for the generalized Chaplygin gas. Rewriting the third 
equation of (6) 



rdS 



(1 + n)- 

l + a 



-d 



;i + H)i+" + 



A 



3 



It then follows to within a factor 5" = (1 + 11) + and r 

Po 

of state Pc — takes the form 



(1+n)- 

l+a 



-A 



A 



l+g 



-. The equation 



(12) 



If we consider that the four velocity has only one component in this model = {N, 0, 0, 0) 
then with the normalization condition we find that 



A 



= /*{ 



Na-^A 



l + a 



1 L i^ + os) 



(13) 



As /i > so (e + 6S) > 0. If we try the canonical methods [24] to find the Hamiltonian 
for this action we will end up with 



(14) 



where = -g^ and Cc is the expression inside second bracket of equation (13). Let us 
now write the super Hamiltonian for the minisupcr space of this FRW model and that is 



1_P 

"24 



a +4F(0)a3^^ + l'^^^ + ) 



(15) 



Analytical solution of the Schrodinger- Wheeler- DeWitt equation formed from this super 
Hamiltonian is very hard to find. So some meaningful approximation can be made for the 
last term. For the early universe scenario we can assume Spl'^"' ^ AaP^^^""^ and using 
this [11] we can expand the expression of the last term of equation (15). 



1 Aa3("+i) 
1 + ^ — + 



-a 



A' 



l + a Spl+^ 2(1 + a)' 



a 



,6(a+l) 



+ . 



Up to the leading order, the super-Hamiltonian takes the form 

1 



,Pl + ST^p, 



24 a AF{(P)a?' 
The canonical transformation 

T = -(1 + a)p~^S^ps , Pt = S'^pe 
would further allow us to simplify the super Hamiltonian and finally 



'H = N 



1 Pi 



24 a AF{(t))a 



:Pl + Pt 



(16) 
(17) 

(18) 

(19) 
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Here Pt is the only canonical variable associated with matter. Now if we can relate any of 
the canonical variables with time, we can have a well defined Hilbert space structure. Let 
us now study the classical dynamics governed by the super Hamiltonian. The Hamilton's 
equation of motion are 



Pa = {Pa,n}=N 



±Pl -J_ 3 
24 a2 ^ 4F{4>)a'^ 



P2 



= 



N 



(20) 



f = {T, n} = N, 

PT = {PT,n} = o, 

with the constraint "H = 0. Actually N acts as a Lagrange multiplier for the system. If 
we see the last two equations of (20) we can see that if TV = 1, then T = t. So T may 
play the role of cosmic time if we choose the gauge = 1. Also Pt = CT(a constant) 
from the last equation. So now with N — 1 and Pt — Ct the third and fourth equation 
of (20) gives 

0, (21) 



0^ 1F'(0) 



2F(0) 



+ 3- 



which on integration yields 



2^(0) = Ca 



where C is an integration constant. The constraint % 
canonical co-ordinates with the help of equation (20) 



(22) 

can be written in terms of 



-Qaa? + ^^F{(l))a? + C^. = 0. (23) 
With equation (22) and equation (23) we can solve for the scale factor and it is 



ait)^ ^(t + CO^-^J% (24) 

J- to satisfy a{t — 



where Ci is an integration constant and we can set it to Ci — y irci^ 
0) = 0. We will now consider the case F{(j)) = 6A^0™, where A and m are constants and 
m 7^ —2 (actually m > give us good results for all values of m). With this and using 
equation (22) and equation (24) we get 



(m + 2) 
6A 



In 



Cat 
t + 2Ci 



2 

m+2 



(25) 
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Figure 1: Left:The classical behaviour of (/) at late times. RightiThe classical behaviour of </> as a 
function of the scale factor. The plots are with numerical values m = 2, A — 0.33, C = 0.375, Ct = 
1, Ci = 1, C2 = 5, C3 = 2. With other admissible values the same nature of the curves are repeated. 

where C2 is a constant of integration. From equation (25) we can see that (j){t =finite) — )■ 
constant. The behaviour of (p at the beginning is not well understood as shown in figure 
1 (left). Let us now find the expression of as a function of the scale factor. With the 
help of equation (22), equation (23) can be written as 



Qaa^ = Ca'-^ + Ct, 
Now dividing equation (22) by equation (26) we get 



F{(j))d(j) 
da 



Ca- 



(26) 



(27) 



Taking only the positive part and integrating we find 

(f){a) = 



3A 




2 

m+2 



+ C3 



(28) 



where C3 is an integration constant. The field blows up at a = and takes a finite 
constant (the integration constant) value for large a. The blowing up of the field is the 
big bang singularity. The classical behaviour of the field as a function of the scale factor 
is shown in figure 1. 



2.1 The Wheeler-DeWitt Equation 

Using the usual quantization procedure let us form the Schrodinger- Wheeler-DeWitt equa- 
tion (also known as Wheeler-DeWitt equation) for our super-Hamiltonian believing that 
the super-Hamiltonian operator annihilates the wave function. So with Pa — ?■ ~ida, 
— 7- —idff), Pt — )■ —idx, h = 1 and 
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we get from equation (19) 



1 



n^{a,(f),t) = 0, 



1 (9^ 
— i- 



24a 4F(0)a3 ^ 



(29) 



(30) 



In order the Hamiltonian operator l-t to be self-adjoint the inner product of any two wave 
functions and \I/2 must take the form [26, 27] 



(*i,*2) = / a'^l^2dad(f), 

J{a,4>) 

and the restrictive boundary conditions being 

9*(a,(/),t) 



or 



^(O,0,t)=O or 



*(a,0,t) = or 



da 

a^(a,0,t) 



= , 



0=0 



(90 



= 0. 



</>=o 



Now we have to apply the method of separation of variables. Putting 

^'(a,0,i)=e^^*77(a,0), 

we get 

where E is the separation constant. If we put 



we will get 



and 



r/(a,0) = x{a)y{(t)), 



+ 



24 



+ 24Ea-^a; + —x = 



da^ ■ " " " ■ 4 

where | is the separation constant. Putting the same expression for F(0) as before 

F(0) - GAV"" , m ^ -2 , 
we get the solution of equation of equation (37) as 



2/(0) = k^+^J(t) 



Ci J_ 



'+2 V m + 2 



m+2 m+2 

(n 2 + c? J -1 I (h 2 



■'+2 \m-\-2 



(31) 

(32) 
(33) 

(34) 
(35) 

(36) 
(37) 

(38) 

(39) 
(40) 
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The solution of equation (38) is also known and is 



x{a) = -\/a 



C3J. Q\/2lBai) + (^^V24eJ 



(41) 



where u — ^'^ and C(i^2,3,4) are integration constants. If we consider C2 = C4 — the 
boundary conditions are well satisfied and the wave function of equation (29) is 

^/^(a,(/),t) = cs-^^fc^ f^\/24Eai) f -^0"^ ) e^-^* , (42) 

^ V3 / "1+2 ym + 2 J 

where C5 is the combination of preceding integration constants. If we want to construct 
the wave packet as a superposition of the eigenfunctions then then we get 

^wp{aA,t) = ^ ^J^_^A{q)B{k)J, [qJ) J_j_ e'^''"'dqdk , (43) 

correct upto a constant factor. Here q — |-\/24E' and A{q) and B{k) are suitable functions 
constructed from the integration constants so that the integration can be carried out 
swiftly. The factor k">^ of equation (42) is absorbed inside B{k). If we choose 

A{q) ^ q''+'^ 6-^"^ (44) 
the integration on q is known [29] and we get 

^wv^\ a4'-, — r / B(k) -^^J 1 2 \ dk , (45) 

where s — ^y — i-^t. At this stage if we approximate y/1 — k'^ ~ 1, then the general pattern 
of the resulting wave packet will remain same [18]. Now if we choose 

B{k) = k^^^^Vl^ , (46) 
the integration can be done [30] and in final form 

, -2m-3 e~47 f ^ <ni+2.\ / tn\ 

= 00 2 — -T^M — 2 47 

(2s)3 -"+2 \m + 2 y 

correct upto a constant factor. In the next section using this wave packet we will follow 
two different approaches to recover the classical solutions of the dynamical variables. In 
figure 2 we have plotted the square of the wave packet against the scale factor and the field 
(f) for different values of time. There we see that the wave packet has two distinct peaks 
and both the peaks have the same height. This can be interpreted as that our universe 
might have evolved from one of these peaks and tunnelled with time. Looking at the figure 
we can also see that the peak height is getting reduced with time but progresses along 
increasing scale factor with a finite spread. Throughout this evolution we also notice that 
the field almost remained constant which is expected. As we have already approximated 
the super Hamiltonian to study the dynamics at a very early universe, so we can say that 
the field approached a nearly constant value within a very short time after big bang. 
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\llf\Ut = 0) |,/r|2(t = 10) 




Figure 2: The behaviour of \'i>wp{a,4',t)\^ at four different times 0,10,30,60 is shown. The plots are 
with numerical values to = 2,A = 4,7 = 0.09375. With other admissible values of these parameters the 
same nature of the curves are repeated. 
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2.2 Expectation values from the Many- World Interpretation of 
Quantum Mechanics 

Using the many-world interpretation of quantum mechanics [28] we can calculate the 
expectation value of the scale factor a and the field 0. Using equation (31) and equation 
(47) one can calculate the expectation value of the scale factor 



< a> {t) 



Straightforward calculation gives 



(48) 



<a>{t)- ^3 



r(|) 



27 + 
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5127 



(49) 



Similarly we can calculate the expectation value for the field (j). In this case 



< > (i) = 



(50) 



which on integration yields 

roo 

J4>=-oo ' 



<<P>{t) 



J 2m+6 , , „ 

\m+2 



dcf) 



J 2m+5 . , „ 



constant. 



(51) 



Prom equation (49) we can see that the expectation value of the scale factor do not 
vanish at t = 0. So it represents a bouncing non-singular universe and may behave like 
a dust dominated classical universe at later times. From equation (51) we see that the 
expectation value of is independent of time. It agrees well with the interpretation from 
equation (25) and with equation (28) (which says that catches a constant value for large 
values of time and scale factor) and the behaviour of as shown in figure 2. 



2.3 Behaviour of the scale factor and the field from the Onto- 
logical Interpretation of quantum mechanics 

Let us now evaluate the expression for the scale factor a and the field using ontological 
interpretation of quantum mechanics. In this formalism the wave function ^ is written 
in the form 

^ = Vle'^ (52) 
and the conjugate momenta corresponding to the variable q is given by 



(53) 
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Portraying the wave function of equation (47) in the form of equation (52) we find 



Bla, (p,t) — - arctan 

^ ^ 3 



3aH 



12872 + U 



9+2 



So with equation (53) 



Pa 



9aH 



P4> 







(54) 



(55) 



12872 + 

Prom equation (9) with our previously chosen gauge N — 1 and with equation (55) we 
get 



a{t) = Oq 



512 o 3 

— 7^ + -t' 
3 ' 2 



and 



0(t) = constant 



(56) 



(57) 



ao is an integration constant. Here also the scale factor is not zero for t = 0, so repre- 
sents a non-singular bouncing universe and 4> is independent of time. So the many-world 
interpretation and the ontological interpretation give us the same behaviour for the scale 
factor and the field. But in general both these interpretations may not predict the same 
behaviour of the variables(seen in anisotropic cosmological models [22]). 



3 Conclusions 

Here we have studied a fiat minisuperspace FRW cosmological model with Chaplygin gas 
and a scalar field coupled to the metric. Schutz's mechanism has allowed us to obtain 
the Schrodinger- Wheeler- De Witt (SWDW) equation for this minisuperspace in our early 
universe. The Chaplygin gas gave us the only remaining matter degree of freedom which 
finally played the role of cosmic time. We have constructed a well behaved wave packet 
from the linear superposition of the wave function of the SWDW equation. While solving 
the SWDW equation we did not bother about the factor ordering of the position and 
momentum operators present in the equation and it is seen that the behaviour of the 
constructed wave packet remains same for other factor orderings. Most interestingly 
we observe two peaks in the wave function, which could be interpreted as, that our 
universe might have evolved from one of those pcaks(analogous to states in quantum 
mechanics). But the choice completely depends upon the initial conditions imposed on 
the field. Here we have observed that the wave packet propagates in the direction of the 
scale factor with time but the value of the field remained constant. We have calculated 
the time evolution of the expectation value of the scale factor and the field using the 
many-world interpretation of quantum mechanics and found that a bouncing and non- 
singular universe. The divergence of the field is also seen at big bang which assumes a 
constant value at later times. We have also calculated the time evolution of the scale 
factor and the field using ontological interpretation of quantum mechanics and found the 
same results. 
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